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ABSTRACT 

Using data of nearby galaxies from the Sloan Digital Sky Survey we investigate whether stellar mass 
(M s tar), central velocity dispersion (er), surface mass density (E), or the Sersic n parameter is best 
correlated with a galaxy's rest-frame color. Specifically, we determine how the mean color of galaxies 
varies with one parameter when another is fixed. When M sta r is fixed we see that strong trends remain 
with all other parameters, whereas residual trends are weaker when E, n, or a are fixed. Overall er 
is the best indicator of a galaxy's typical color, showing the largest residual color dependence when 
any of the other three parameters are fixed, and M star is the poorest. Other studies have indicated 
that both the halo and black hole properties are better correlated with a than with M sta r, E or n. 
Therefore, our results are consistent with a picture where a galaxy's star formation history and present 
star formation rate are determined to some significant degree by the current properties and assembly 
history of its dark matter halo and/or the feedback from its central super massive black hole. 
Subject headings: galaxies: fundamental parameters — galaxies: formation — galaxies: kinematics 
and dynamics — galaxies: statistics 



1. INTRODUCTION 



established that the stellar populations 
in the nearby Universe correlate with 



It is well 
of galaxies 

their luminosities and masses, such that the stars in 
more massive g alaxies are on average older and more 
metal-rich fe.g.JBower et al.lll992t IBlanton et al.ll2003at 
IKauffmann et all l2003f T~ The old ages of the most lu- 
minous galaxies are somewhat puzzling, as it implies 
that there must be some mechan ism responsible for shut 
ting off their star formation fe.g..lKauffmann &: Haehnel 
2000; iCroton et aLll2006t iNaab et al.ll2007t iKeres et al 
20051 ). Another puzzle is that the correlation between 
star formation history and mass is complex: low mass 
galaxies are typically younger and high mass galaxies are 
typically older, with a bimodal tr ansition region at ~ 3 x 
10 10 M Q (|Kauffmann et al J 120031) . This bimodality and 
the associated transition mass scale have been the sub- 
ject o f intense debate i n the past years (e.g.jBundv et al.l 
12001 |Pen7eTaT]l20Tot iBrammer et "aLll2011h . 

An intriguing possibility is that luminosity and mass 
are not the "right" parameters to interpret galaxy evolu- 
tion, and that there exist other parameters that show 
more straightforward correlations wit h stellar popula- 
tion parame ters. As demonstrated by IKauffmann et all 
(|2003ll2006T ), the star formation histories of galaxies show 
less scatter when the structure of galaxies is taken into 
account. In particular, the ages and star formation rates 
of galaxies are better correlated with their surface den- 
sities E (which is proportional to M.*jr\ with r e the 
size o f the galaxy) than with mass alone. iFranx et al.l 
(2008) find that the strong correlation between color and 
M*/rg (and M*/r e ) persists all the way to z ~ 2. Sim- 
ilar tre nds have bee n found for v elocity dispersion and 
for the lSersid (|1968h index n (e.g.. IBlanton et al.ll2003al : 
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Although it seems clear that the correlation between 
color and mass is weaker than the correlations between 
color and E, a, or n, it is not clear which of these pa- 
rameters is the best predictor of a galaxy's color. This 
is important to establish as it provides information on 
the physical processes that govern galaxy evolution. As 
an example, if n b est correlates with color (as suggested 
bv lBell et atll201lD it may imply that the merger history 
determines the star formation history, whereas if a is the 
key parameter it would suggest that the star formation 
history is influenced (or even determined) by the dark 
matter halo or the central supermassive black hole. 

In this Letter we investigate which of the four parame- 
ters M star , E, cr, or n is most closely correlated with the 
star formation history, as parameterized by the color. 
We determine this by fixing each parameter in turn and 
measuring to what extent the color depends on the other 
three parameters. The homogeneous, large datasets re- 
quired for this study are now available from the 7th 
data release of the Slo an Digital Sky Survey (SDSS DR7; 
lAbazaiian et al.H2009fh 

2. DATA 

The galaxy data used in this analysis are gathered from 
the s eventh data release of the SDSS (|Abazaiian et al.l 
2009). We begin with the Large Scale Structure sam- 
ples of the DR7 NYU Valu e Added Galaxy catalogue 
(VAGC Bla nton et alll2005l) . The sample we use has 
an r band magnitude range of 14.5 < r < 17.6. In 
addition the NYU V AGC gives k-correct ed (to z=0.1) 
absolute magnitudes (Blanton et al. 2003a), velocity dis- 
persion measurements from the Princeton Spectroscopic 
pipeline, and circularized sersic fits for each galaxy all of 
which we make use of in this analysis. 

For estimates of the stellar mass we make use of the 
MPA/JHU DR7 value added catalog which provides stel- 
lar mass estimates based on stellar population fits to the 
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Figure 1. The u — r color distribution as a function of M sta r, 
a 2 , E and n 2 , for SDSS DR7 galaxies with 0.04 < z < 0.113. 
The u — r color is positively correlated with all four parame- 
ters. The correlation with M st ar has large scatter, particularly 
near M star ~ 3 x 10 10 Mq. The dotted lines show the cuts we 
apply in M star and a to define our primary sample. 

SPSS photometry (jKauffmann et al.l l2003t iSalim et al.l 
[2001 . The overlap between the MPA/JHU and NYU 
VAGs is close to but to quite 100% and so we remove 
the regions where they do not match from the analysis. 
We also remove any region of the survey that has a spec- 
troscopic completeness less than 70%. This leaves an 
area of 7640 deg 2 and a total sample of 521,313 galaxies. 

The SDSS velocity dispersions are measured within 
the 3" diameter SDSS fiber. We correct to a com- 
mon aperture of one eighth of an effective radius (r e ), 
the central velocity dispersion, using t he relation op = 
o~av(& r av /r e ) 066 where r ap = l."5 (jCappellari et al.l 
2006) . r e is taken from the best fitting circularized Sersic 
profile fit. 

Throughout we use u — r color from the if-corrected 
NYU VAGC absolute magnitudes. As already stated 
these are corrected to z = 0.1; although they are not 
quite it - r at rest we will refer to them as u — r col- 
ors. We also make use of the morphologic al classifica- 
tions from Galaxy Zoo (|Lintott et al.ll201l[) which pro- 
vides multiple visual classifications for each galaxy in the 
SDSS spectroscopic sample. The parameter we use is the 
probability that a galaxy is an edge on disk (P e dge) (see 
iLintott et aT]|201lL for details). 

3. DEPENDENCE OF COLOR ON M star , <r, E, AND N 

We wish to understand how the typical color of a 
galaxy depends on its stellar mass, central velocity dis- 
persion, surface mass density, and Sersic profile. We be- 
gin in Figure [1] were we show the u — r color distribution 
for all SDSS galaxies with 0.02 < z < 0.11 as a func- 



tion of each parameter. We have applied to this figure 
(and all subsequent figures) a V/Vmax weight for each 
galaxy to correct for the varying stellar mass complete- 
ness limit with redshift. We plot the color distributions 
as a function of log(M star ), log(S), log(cr 2 ) and log(n 2 ). 

We have chosen these units as they are all approxi- 
mately linearly proportional to one an other. This is 
illustrated in Figure [2] were we show the relationships 
between these four parameters. The contours show the 
full distributions and the points show the mean of y in 
bins of x (black) and y in bins of x (red). Whilst we 
may expect £ and a 2 to be approximately linearly pro- 
portional to M s tar and each other, we find that n 2 shows 
the same trends over a broad range in those parameters. 
Therefore we can meaningfully compare the predictive 
power of n 2 to that of the other three parameters. 

Figure [T] shows that there is a clear color dependence 
on each parameter but with significant scatter. It also 
illustrates, as has been previously found (e.g., Kauffmann 
ct al. 2003), that M star is a relatively poor discriminator 
of a galaxy's color, particularly around the "transition 
mass" of ~ 3 x 10 10 Mq. It is important to note that 
distributions become noisy at both low M star and low 
a. The number of galaxies becomes increasing small at 
low masses, due to the apparent magnitude limit of the 
SDSS, which also strongly affects the sample at low a 
due to the tight correlation between the two parameters. 
a measurements also become more uncertain below 65 
km/s due to the resolution of the SDSS spectra. For 
these reasons we limit the remainder of the analysis to 
galaxies with M sta r > 10 10 M Q and a > 65 km/s. 

We also remove edge on disk galaxies from our sample 
to m inimize the influe nce of dust on our measurements 
(e.g. JPatel et al.ll2011l) . although we expect the dust con- 
tribution to be low in the local universe. We make a cut 
using the Galaxy Zoo P e dge parameter limiting it to be 
< 0.3. We note that this cut changes our results very 
little and has no effect on any of our conclusions. 

4. WHICH PARAMETER CORRELATES BEST WITH 
COLOR? 

In order to determine which parameter shows the 
strongest correlation with color we determine whether 
there remains any color dependence on M s t ar , S, er, and n 
when each of these parameters are held fixed. We divide 
our parent sample ( 0.02 < z < 0.11, M star > 10 10 M Q , a 
> 65 km/s, Pedge < 0.3) into a series of samples selected 
to have narrow ranges in each parameter. We select bins 
of 0.05 in log(M star ) and log(S), 0.025 in log(er) and 0.2 
in n. We then calculate the mean u — r color of galaxies 
in each of these narrow binned samples as of function of 
the other three parameters, where the mean is calculated 
in bins of 300 galaxies. 

4.1. Residual Correlations 

We show the resulting relationships between mean 
color and M star , S, n and a in Figures [3] and |4j Each row 
shows a pair of parameters with the binning and abscissa 
parameter switched. That is, in each row the left panel 
shows the effect of varying parameter 1 at fixed param- 
eter 2, and the right panel shows the effect of varying 
parameter 2 at fixed parameter 1. The left and right 
panel therefore essentially show the same information, 
but highly the trends in a complimentary way. If the 
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Figure 2. The relationships between M sta r, S, <r 2 and n 2 for all SDSS DR7 galaxies with 0.02 < z < 0.11. Galaxies are V/Vmax 
weighted to correct for the redshift dependent stellar mass completeness limit. The black points show the mean of y in bins of 
x containing 500 galaxies, where as the red points show the mean of x in bins of y. We can see that over a significant range of 
each parameter they are all approximately linearly proportional to one another. 



trends of all of the binned samples lie on top of each 
other in any of the plots it would mean that the color 
would be independent of the parameter used for the bin- 
ning. Similarly, if the abscissa parameter is unimportant 
each of the individual binned trends will be flat. 

Turning to the first pair of parameters at the top of 
Figure |31 £ and M star , we can see something close to 
this extreme situation. When M sta r is held fixed (left 
panel) there remains a clear trend with S, such that 
higher £ galaxies are redder. The trend is particularly 
strong at low M sta r and gradually decreases at the high- 
est masses, resulting in some spread but a generally low 
dispersion between the M star bins. The same trend is 
visible when £ is held fixed; there is only a very weak 
dependence of mean u — r color on M star with more mas- 
sive galaxies being redder. The individual £ bins are 
well separated, showing the strong color dependence on 
£, and are generally parallel. It appears that S is a 
better indicator of a galaxy's color than its stellar mass 
is. 

We next consider M star and n, shown in the middle 
panel of Figure G3 For low n (< 2.5) there is a strong de- 
pendence of the color on n at fixed M star , but for higher 
n the trends with M star are all the same with no n depen- 
dence. At all n there is a trend for more massive galaxies 
to be redder on average, although this is very weak for 
galaxies with log(M star ) < 11. Again M star appears a 
poorer indicator of a galaxies color than n, although this 
is less true at high n or high M star . 



The bottom panels of Figure[3]concern £ and n. At low 
n there is a very strong relationship between color and 
£ with higher E galaxies being redder. Whilst this trend 
remains for all n the slope of this correlation reduces as 
n increases. There is a similar dependence of color on n 
at fixed E, although the relation becomes very weak at 
the highest S. 

In Figure [4] we show how the color depends on a at 
fixed M s t ar , E, and n on the left and how color depends 
on M star , E, and n at fixed a on the right. It is striking 
how tight all of the trends are in the left panels and how 
separated they are in the right panels. This indicates 
that the mean galaxy color is more strongly dependent 
on a than M sta r, E or n. This is particularly true at 
a (> 200 km/s) where all the trends lie almost completely 
on top of each other on the left and are very close to 
flat on the right. At lower a interesting trends emerge; 
perhaps surprisingly the mean color becomes bluer as 
M s tar increases, with this trend increasing as a decreases. 
This probably reflects an increasing disk component in 
more massive galaxies at fixed a. At low a (<175 km/s) 
a trend emerges with E such that galaxies with higher 
E at fixed a have redder colors, with the trend becoming 
more significant as a decreases. There is also a trend 
with n when n < 2.5 such that galaxies with higher n 
are redder. There is no color dependence on n at fixed 
a for higher n galaxies. Whilst some trends with M s tar, 
E and n emerge in some regions of the parameter space 
there is always a strong dependence of the mean color on 
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Figure 3. The relationship between mean u — r color and mass surface density at fixed stellar mass (top left), stellar mass at 
fixed mass surface density (top right), stellar mass at fixed sersic n (middle left), sersic n at fixed stellar mass (middle right), 
mass surface density at fixed sersic n (bottom left) and sersic n at fixed mass surface density (bottom right). 
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Figure 4. The relationship between mean u — r color and velocity dispersion at fixed stellar mass (top left), stellar mass at 
fixed velocity dispersion (top right), velocity dispersion at fixed mass surface density (middle left) and mass surface density at 
fixed velocity dispersion (middle right), velocity dispersion at fixed sersic n (bottom left) and sersic at fixed velocity dispersion 
(bottom right). 
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a over the full range of the other parameters. 

4.2. Quantifying the Residual Correlations 

We show in Figure [5] an attempt to both simplify and 
quantify the trends that are displayed in Figures |3] and 
[4j In narrow bins of one parameter we calculate the dif- 
ference in the mean u — r color of the galaxies lying in 
the lowest and highest 10 percentiles of the other three 
parameters. Since the size of the range of the second 
parameter may vary with the first and between the dif- 
ferent parameters we calculate this color difference per 
unit log(M s tar), log(S), log(er 2 ) and \og{n 2 ). For exam- 
ple the blue points in the top left panel of [5] show the 
difference between the mean u — r color of galaxies with 
the 10% highest E and 10% lowest E per unit log(E) in 
narrow bins of a. This is essentially the vertical scatter 
in the center left panel of Figure 0] or the gradient of the 
individual trends in the center right panel. 

The color gradients shown in Figure [5] reinforce the 
conclusions we have already drawn. When a is fixed 
(top left panel) the magnitude of the color gradient is 
always less than 0.5, showing that there is only a weak 
dependence of color on M sta r, E or n at constant a. This 
is particularly true at high a (> 200 km/s) where the 
gradient is essentially zero for both E and n. When the 
other three parameters are held fixed the color gradient 
with a is always the largest and is greater than 0.5 re- 
gardless of the values of the other parameters. Clearly 
there is a significantly larger dependence of the color of 
galaxies on a than on M sta r, E or n. 

Figure [5] again confirms the negative color gradients 
with M s tar at low cr (red points in top left panel), such 
that more massive galaxies are bluer when a is fixed. We 
also see negative color gradients with n for the highest 
E galaxies. Otherwise it is always the case that galaxies 
with higher M sta r or n are redder, and always the case 
that galaxies with higher cr or E are redder when the 
other parameters are fixed. 

5. DISCUSSION 

The central result of this paper is that the colors of 
galaxies depend more strongly on a than on M sta r, E, or 
n. We have demonstrated this by examining, for each of 
these parameters, how strong residual correlations with 
the other three parameters are when the parameter under 
consideration is held fixed. At fixed a residual trends 
with other parameters are weak, and when any of the 
other parameters is fixed there are strong residual trends 
with a. 

Recently iBell et ahl (|2011[ ) studied how the fraction of 
passive (i.e., not star-forming) galaxies depends on a sim- 
ilar set of parameters and reached the conclusion that 
n was the best indicator of the passive fraction. This 
seems to be at odds with our findings where a shows the 
strongest color trends. Th e difference may stem from the 
fact that IBell et ail (|20ll use M star /Re to approximate 
velocity dispersion. Whilst there is a correlation between 
cr and M s tar/-Re, which is improved if corrections based 
on n are applied (jTavlor et alj|2010h . there is significant 
scatter (~0.1 dex). This may be the cause of our differ- 
ing conclusions and illustrate the importance of directly 
measured dynamical properties. 

Ou r r esults, and t h ose of ot hers Ce.g..lKauffmann et all 
[200l iFranx et all [20081 Ivan Dokkum et all l201lT 



IBell et all 120111) . call into questi on whether "mass 
quenching" (e.g., iPeng et all I2010D and other mass- 
driven effects are actually manifestations of underlying 
trends with velocity dispersion. The velocity dispersion 
may in turn reflect a yet more fundamental parameter. 
It is known to correlate well with central black hole 
mass (e.g., iFerrarese fc Merrittl l2000t IGebhardt et all 
120001) and also wit h the properties of dark matter halos 
(Wak e et all 120121 Taken together, all these results 
suggest a simple scenario where the mass of the dark 
matter halo determines both the central black hole mass 
and the star formation history of galaxies. 

This study can be extended and improved in many 
ways. The dispersions are currently corrected to a com- 
mon aperture with reasonable assumptions, but it would 
be very useful to measure radial trends in dispersion in a 
systematic way. This is particularly relevant for low mass 
galaxies and star forming galaxies, as they have signif- 
icant disks which presumably dominate at large radii. 
Blue disks may well be the cause of the peculiar fact 
that, at low cr, more massive galaxies are bluer at fixed 
dispersion (see Fig. 4) . Modeling of the effects of various 
physical processes (such as merging) on the velocity dis- 
persion may help us understand why velocity dispersion 
is so well correlated with many aspects of galaxies. One 
possibility is that cr may well be indicating both the halo 
mass (or other halo prope rties) , in a similar or more pre- 
cise manner than Mstar (|Wake et al.M2012f) . and at the 
same time be an indicator of the relative bulge to disk 
components. So at fixed M s tar a higher a galaxy has a 
larger bulge to disk ratio and so is redder, whereas at 
fixed n a higher a galaxy typically occupies a more mas- 
sive dark matter halo and so is also redder, a is then 
the best of the four parameters at encapsulating both 
color dependences (halo mass and bulge to disk ratio) as 
illustrated by the tightness of the a - M s tar and cr - n 
relations shown in Figure [5J 

Finally, it would be very interesting to do similar sys- 
tematic studie s at earlier co s mic ep ochs, extending those 
pionee red by IFranx et all (|2008f) and IBezanson et all 
(pOTlT) . which will give further insight into the physical 
parameters which drive galaxy formation and evolution. 
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